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Abstract 

It has been recently found that a number of systems displaying crackling noise also show a remarkable 
behavior regarding the temporal occurrence of successive events versus their size: a scaling law for the 
probability distributions of waiting times as a function of a minimum size is fulfilled, signaling the existence 
on those systems of self-similarity in time-size. This property is also present in some non-crackling systems. 
Here, the uncommon character of the scaling law is illustrated with simple marked renewal processes, built 
by definition with no correlations. Whereas processes with a finite mean waiting time do not fulfill a scaling 
law in general and tend towards a Poisson process in the limit of very high sizes, processes without a finite 
mean tend to another class of distributions, characterized by double power-law waiting-time densities. This 
is somehow reminiscent of the generalized central limit theorem. A model with short-range correlations is 
not able to escape from the attraction of those limit distributions. A discussion on open problems in the 
modeling of these properties is provided. 
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I. INTRODUCTION 



In the words of Sethna et al, "crackling noise arises when a system responds to changing 
external conditions through discrete, impulsive events spanning a broad range of sizes" yj]. Oper- 
ationally, a broad range of sizes essentially means that the size s of the events fluctuates following 
a power-law distribution, D{s) \/s l+ $ , where D(s) is the probability density of s and 1 + /3 is 
the exponent. This is quite remarkable, as power laws signal the absence of characteristic scales, 
in this case of event sizes Q2Q . It is also implicit that the driving that makes the external conditions 
change is relatively small and smooth, quite different from the resulting bursty response; therefore, 
crackling noise signals a highly nonlinear behavior. 

Although these ideas have been developed within the physics of condensed matter yj], natural 
hazards show perhaps the largest number and best illustrations of crackling noise yfl, including 
J^Q], landslides and rock avalanches |3, 6], volcanic eruptions [Q], rainfall Jgl, hur- 
solar flares [Jil 11], the activity of the magnetosphere [12], and perhaps meteorite 



earthquakes 
ncanes [91, 



impacts [13] (provided that the Earth moves "slowly" through space). In other catastrophic phe- 



nomena, as forest fires 1114 , 



15|] or the extinctions of biological species II 1611 . the events span a 



broad range of sizes but it is not clear if or in which cases they are power-law distributed. 

Beyond the geosciences, notable examples of crackling noise arise in physiology and human 
affairs, like neuronal firings [fTvll . epileptic seizures Jisll and appearances of words in texts or 
speech [ JjJ] (if the "size of a word" is measured by its rarity, i.e., its position in a ranking of 
frequencies). Note that most examples of crackling noise arise in systems with a high degree of 
complexity, characterized by an enormous number of degrees of freedom that interact between 
them. 

Crackling noise can be considered as the most important and apparent property of systems dis- 
playing self-organized criticality (SOC). This concept goes one step beyond, and proposes that the 
origin of the power-law distribution of sizes (i.e., the hallmark of crackling noise) is the existence 
of a critical point (analogous to those of equilibrium phase transitions, which are well-known to 
have scale-invariant properties) to which the dynamics of the system is attracted by means of a 
feedback mechanism that balances driving and dissipation [2(], 21]. The paradigmatic example is 
a sandpile over an open support to which grains are slowly added: when there are few grains, the 
pile is flat and grain dissipation at the border is low, then the pile grows; in contrast, when there are 
too many grains they easily travel through the system and border dissipation is large, so the slope 



of the pile decreases. At the end (in the attractor), the slope fluctuates around a critical value that 
balances the input and output of grains, and this state should have scale-invariant properties, i.e., 



power-law statistics, as in equilibrium critical phenomena. The behavior of real-wor 



d sandpiles is 



22b. Then, SOC 



more diverse than what the SOC picture suggests, but that is a different story 
is a plausible physical mechanism for the emergence of crackling noise, but this does not preclude 
that other mechanisms could also lead to crackling noise [23]. 

In practice, although it is very simple to test if a system displays crackling noise (just measuring 
the size of the events and calculating their distribution, checking that the driving on the system is 
slow and smooth), it is not so easy to demonstrate the existence of SOC, as one would need to 
measure the relations between the internal variables of the system and their fluctuations, and they 



should behave in the same way as the equivalent ones in an equilibrium phase transition [24] 



H. SCALING LAW FOR WAITING-TIME STATISTICS 

Neither the above definition of crackling noise nor the usual studies of SOC pay too much 
attention on how the discrete and impulsive response events that the system develops occur in 
time. Is there a unique dynamical process that defines these behaviors? And in that case, is it 
periodic? Is it chaotic? Is it random? How does the dynamics reflect the scale-invariant properties 



of such systems? |25|]. 



The case of earthquakes exemplifies our poor understanding of the dynamics of this kind of 
processes [25, 2(3]. On the one hand, there is a widespread belief in the notion of characteristic 



earthquakes: the strongest events that a single fault segment is able to generate are always almost 
the same (same epicenter, same size, same focal mechanism) and should occur at regular times 



[27, 



28fl. On the other hand, for extended regions, it is often assumed that mainshocks come in 



a total random way, i.e., from a Poisson process, and aftershocks follow a different process [29] 



Kagan has strongly argued against these simplistic views, showing evidence of time clustering in 



3 



31b 



earthquake occurrence, just the behavior opposite to the characteristic-earthquake concept [3 

Recently, it has been found that some of the systems mentioned above as prototypical of SOC 
or crackling noise display a remarkable temporal behavior. For such systems, let us consider the 
waiting time, also called recurrence time, inter-event time or inter-occurrence time; this is the 
time between consecutive events above a size threshold. So, we take into account only events 
whose size s verifies s > c, where c is the threshold value (but notice that for instance in the case 



of earthquakes one does not distinguish between foreshocks, mainshocks, and aftershocks). This 
defines a set of occurrence times, tf, denoting the occurrence of the z-th event above c, from i = 
to N c . As each event is characterized by a unique occurrence time it is assumed that the duration 
of the event is very short in the scale of observation, and therefore the process can be described 
as a stochastic point process |32, 33|. In addition, each event is also characterized by its size, so 
the process can be considered a marked point process, the size being a "mark" added to the time 



occurrence (we do not consider spatial degrees of freedom in this paper, but see Refs. 1134 



any case, the waiting times for events with s > c are obtained straightforwardly as xf = tf — tf 



350). In 

c 

i-V 



with/ = 1,2... N c . 

I — I 

The key element of analysis was introduced by Bak et al. [|36fl , by means of a systematic study 
of the statistics of xf as a function of the size threshold c. Although the rise of c only eliminates 
some values of the occurrence times, leaving the rest unaltered (i.e., tf — > tj , with j < i and c < c'), 
the waiting times are changed in a more complicated way as they add in a variable number to give 
rise to the larger (or not) new waiting times. Usually, for the type of systems that have been studied 
so far, the waiting times show a large variability, and the best characterization on these processes 
is by means of the waiting-time probability density 1371]. 

What has been found is that for many such systems, these probability densities verify a scaling 
law. If for events with s>cwe denote the waiting-time probability density by D c (x) and the mean 
waiting time by x c (f c = Jq xD c (x)dx), the scaling law can be written as 

D c (x)=F(x/x c )/x c , (1) 

where F is a scaling function, independent on c. This means that the shape of the distribution 
is independent on the scale given by x c (which obviously is determined by the threshold c); in 
other words, when the waiting time is measured using as a unit its mean value, the results are 
independent on the value of c, which implies the existence of a self-similarity in the process. We 
will argue in the rest of this paper that this is quite a remarkable result in general, difficult to justify 
with the use of simple stochastic models. 

In the case crackling noise or SOC systems, the mean waiting verifies x c <x \j D(s)ds c& , 
if /3 > 0, and substituting in the scaling law, 

*>c(T) = \F (Vj , (2) 



cP \cP. 

where F is the scaling function F incorporating the factor of proportionality between x c and c$ . 
Written in this form, the scaling law turns out to be a particular case of the condition of scale 



invariance for functions with two variables, T and c |2fl. So, although for one variable the condi- 
tion of scale invariance yields a power law (for instance, for s we have D(s) l/s i+ P), for two 
variables, like s and T, scale invariance leads to Eq. ©, with F an undetermined function. 



Crackling-noise or SOC systems showing this behavior include earthquakes 



380, fractures 



solar flares Qll ]_. 



with previous belief 1142 , 



iterary texts [41], or some paradigmatic sandpile models, in contrast 



43|] . But this property is shared by other systems for which its crackling 



440, forest fires LI 5] and 



noise nature is in doubt, as printing requests in a computer network 
tsunamis Ii450 (although the latter seem to be power-law distributed [14611 . certainly they are not 
slowly driven; rather, they are cracklingly driven by undersea earthquakes and landslides). Even 
systems that do not crackle, as diverse climate records (temperature, river levels, etc.) [|470, or 
systems for which the crackling behavior is in the derivative of the response signal, as financial in- 
dices |4s|, verify a scaling law as Eq. © when the threshold is large enough that the events above 
it become extreme events. The corresponding scaling functions come in a variety of functional 
forms; are there any preferred types? 



III. MODELS FOR TIME-SIZE SCALE INVARIANCE OF EVENT OCCURRENCE 



A. Marked Poisson Process 

Which is the meaning and depth of the scaling law (OQ)? Certainly, a marked Poisson process 
trivially fulfills it. This is a marked point process in which the occurrence times follow a Poisson 
process, and the sizes of the events (the "marks") come from a random distribution independently 
on occurrence times and other sizes. Its simulation is very simple, with independent identical 
exponentially distributed waiting times and independent identical power-law distributed sizes (in 
the case of crackling-noise systems). 

Indeed, a Poisson process is completely characterized by its rate, let us say, 10 events per hour. 
If we now raise the size threshold in such a way that half of the events are eliminated and half 
of them survive, this is equivalent to a random thinning of the events with a thinning probability 
equal to 1/2 (in which any event has the same probability of being eliminated, independently of 
the rest), due to the fact that the sizes of the events are uncorrected. So, we end with a Poisson 
process of rate equal to 5 events per hour. 

It is well known how to show this more rigorously. Consider that events are removed from a 



marked Poisson process with a probability q, and are kept with probability p = 1 — q; then, the 
probability that the number of events N' that survive in a time interval of length A is equal to k is, 

oo 

Pr[JV' = k] = £ Pr[N' = k\N = n] Pr[N = n] , 

n=k 

where Pr denotes probability, "|" conditional probability, N is the number of events in the interval 
prior to thinning, and n counts all the possible values of N. By hypothesis, the original process is 
Poisson of rate X, so Pr[N = n] = e~ v v n /n\, with v = A, A, and by the uncorrected nature of the 
process Pr[N' = k\N = n] is given by the binomial distribution, 

Pr[iV' = k\N = n} = 

Substituting both above, 

Fr \N> = k]= e-v p - f -^v" = l^£ e -v f Ml! - M e -v(W) _ e - P A_Pvl 
1 1 k\^ k {n-k)l k\ h( n ~ k V- kl k\ ' 

which defines another Poisson process of rate pX . If we rescale the new rate X' = pX as X' — ► X' /p 

we recover precisely the original Poisson distribution. Note that p is given by p = Pr[s >c'\s>c], 

and in our context it turns out that p = [c/c')^ . 

So, could the trivial marked Poisson process explain the scaling law Eq. (OQ)? Certainly not, 

as none of the known examples mentioned above are characterized by an exponential scaling 

function. We should go beyond this trivial explanation. 




B. Marked Renewal Process 



The shapes of the scaling functions found for the real data mentioned above are rather diverse, 
including the gamma distribution, the stretched exponential, and the power law for large times. It 
seems necessary to incorporate this shape into the point process modeling those systems. The most 
straightforward way to do this is through a renewal process, which is characterized by independent 
identically distributed waiting times, following a specific distribution [32]. If we add to this model 
independent identically distributed sizes we end with a process that we may call marked renewal 
process. Note that there are no correlations whatsoever in this process, but there is a memory of 
the last event if the waiting-time distribution is not exponential. 

The probability distribution of waiting times for events with s > c' can be obtained from the 
probability distribution for those with s > c, if d > c. The idea is the same as in the previous 



subsection but we will use the waiting-time representation rather than the count-number represen- 
tation of the process. The same steps as in Ref. [49] will be followed, although the case here is 
simpler. 

We start using the survivor function, S c >(x) = Pr[ waiting time > x for events with s > c'} = 
f™ D c /(x)dx. If an event of size so > c' has taken place, the next one with s > c' can happen in a 
variety of ways, depending on the number of events with c < s < c' in between. So, we can write, 

oo 

S c ,{x) = £ Pt[t^ > T, Jl< c',. .. Sj-l < c', Sj > c'} = 

£ Pr[T (;) > T 1 51 <</,... Sj-l < c', Sj > c] ■ Pr[5i < c', . . . Sj-l < c', Sj > c] = 

j i 

oo 

Pr[T (;) > x] ■ Pr[5i < c'] ■ ■ -Vx[sj-i < c'] ■ Vr[sj > c'] = 

j 1 

where Pr denotes probability, "|" conditional probability, and the j— th return time is defined, for 
events with s > c, as = ? ( c — tf_j, that is, as the elapsed time between any event and the j— th 
event after it (naturally, the first return time is the waiting time). The conditions on Pr[r^ > t] 
are eliminated because waiting times are independent on sizes. As in the previous subsection 
p = Pr[5 > c' | s > c] = Pr[s > c'] and q = 1 — p = Pr[s < c' \ s > c] = Pr[s < c'] (the condition s > c 
is always implicit, if it is not explicit). Therefore, 

oo 

Sj{x) = Ypq j - x H* U) >A- 

If in this equation we derive with respect x we obtain the probability densities of the return 
times; as the waiting times are considered independent on each other, we use that the j— th-return- 
time distribution is given by j convolutions of the first-return-time distribution (denoted by the 
symbol *) to get 

oo 

DA?) = £p<7 J ' -1 [A:(T)n =pD c (x) + qpD c (x) *D c (x) +q 2 P D c (x) *D c (x) *D c (x) + ■ ■ ■ (3) 

where the exponent * j means that D c (x) is convoluted with itself j times. It is convenient to look 
at Eq. © in Laplace space, where things are simpler, then D(s) = [q e~ ST D(x)dx, which is a 
(moment) generating function, and the convolutions turn into simple products, 

oo 

D c ,(s)=pD c (s)Y,<l J -~ l [Dc(s)} j - 1 = P D c {s)+qpDl(s)+q 2 pD 3 c {s) + --- (4) 

7=1 



As q and D c (s) are smaller than one (this is general for generating functions), the infinite sum can 
be performed, turning out that 

DM = J^L. (5) 
1 -qD c (s) 

Equation © describes the effect of rising the threshold on the waiting-time distribution. The 
next step is the scale transformation, which puts the distributions corresponding to c and c' on the 
same scale. We will obtain this by removing the effect of the decreasing of the rate, which, is 
proportional to p, so, 

D c ,(t) - p-'D^T/p), (6) 

and in Laplace space we get 

Dj(s) - D d {ps). (7) 

Therefore, the combined effect of rising the threshold plus rescaling leads to a transformation T 
that acts on the original distribution, 

TDM = T^^T. (8) 
1 -qD c (ps) 

We are very interested in the fixed points of this transformation, which are obtained by the 
solutions D* (s) of 

TD*(s)=D:.(s) 1 

where * now means fixed point; The previous fixed-point equation is totally equivalent to the 
scaling law (QQ)- Introducing the variable co = ps and substituting p = co/s and q = 1 — co/s in the 
fixed-point equation we get, separating variables, 

_J_-I = ^ 1 = 1. (9 ) 

sD*(s) s coD*((o) a A' 

where we have made both functions equal to an arbitrary constant due to the fact that p and s 
are independent variables and so s and co are; then, the only way in which the equality could be 
fulfilled, for all s and CO, is that the function is a constant 1/A. The solution is then 

d *m=ttW <10) 

which is the Laplace transform of an exponential distribution, 

D* c (x) = Xe- Xx . (11) 



The dependence on c enters by means of X, as A -1 = f c . Note that this demonstration includes the 
one on the previous subsection, showing that the marked Poisson process displays a scaling law 
for the waiting-time distribution, but in this case we have achieved a more general result, as the 
marked Poisson process is the only marked renewal process which can fulfill such a scaling law 
(when the rescaling is done with the mean waiting time f c ). 

We can go one step beyond and demonstrate that the marked Poisson process is an attractor for 
the broad family of marked renewal processes for which the mean waiting time f c is finite. The 
iterative application of transformation T with a finite probability p is equivalent to the limit p — > 

in Eq. ©. Expanding that equation up to first order in p, using D c (ps) = 1 — f c ps H , yields 

pD c {ps) p(l-i c sp-\ ) 



TD c (s) 



l-qD c (ps) l-(l-p)(l-f c sp + - 
1 — t c sp H 1 — T c sp H 1 



+ • 



p-T c sp-\ ) l+T c S-\ 1+T C S 

which indeed corresponds to a Poisson process when p — > 0. 

This result illustrates the strange particularity of the scaling relation (Q~|): among the infinite 
number of probability distributions with a finite mean that can define a marked renewal process, 
only one type, the one with exponentially distributed waiting times, fulfills the scaling law. The 
results can be put in the language of the renormalization group. Indeed, the first part of the pro- 
cess, called thinning, where the threshold is raised from c to c', corresponds to a decimation of 
events. This is analogous to the renormalization of the Ising model, where some portion of the 
spins are eliminated 0, 13]. The second part of the process correspond to a change of scale in 
time, which is equivalent to the change of scale in real space renormalization. A third step, the 
renormalization of the field, is not necessary for the purposes of computing waiting-time statistics. 
So, the transformation T can be considered a renormalization transformation, and we have seen 
how a renewal process (with a finite mean) renormalizes into the trivial Poisson fixed point. So, 
for all the renewal processes of this kind (except for a set of zero measure, given by the Poisson 
process) one expects a change under renormalization, and not scale invariance. This is one of the 
reasons why the existence of the scaling law (Q]) is so intriguing. 



C. Marked Renewal Process without a Finite Mean Waiting Time 



What happens for renewal processes whose waiting-time density does not have a finite mean? 
Obviously, a rate cannot be defined as the inverse of the mean, nevertheless, still it is possible 



to follow an approach that makes sense. The first part of our transformation, in which the size 
threshold is increased from c to c', does not change ©; however, the rescaling with the mean 
cannot be applied. We will use as a rescaling parameter p = Pr[s > c' \ s > c\, but in contrast to the 
previous case we will raise p to some power r in Eqs. © and ©. This is equivalent to seek for a 
scaling law of the form 

D c (T)=R r c F(R r c T), 

where the rate R c is understood as the number of events per unit time in the time window under 
consideration. 

The transformation T [Eq. ©] then becomes 

TDAs) = / D i P ^ v (12) 
1 -qD c (p r s) 

In the same way as before, the fixed point equation leads to 

1111 



> l l r D* c {s) s l / r (0 1 / r D*((o) (O l l r 



= a, (13) 



whose solution is 



1 +as a 



with a = l/r. Only for < a < 1, i.e., r > 1, this function represents a probability distribution; 
this is so because for other values of a the expansion of D* (s) does not correspond to the expansion 
of a generating function. However, when < a < 1 a finite mean does not exist. 

Indeed, for small s, D*.(s) = 1 —as a ; this corresponds, if < a < 1, to the Laplace transform of 
D*(t) = A/z l+a , for large T, with a = —AT(—a) and T(—a) the gamma function of —a [51]. On 
the other hand, the behavior for large s is D*.(s) = 1/ (as a ), and by means of a Tauberian theorem 
the limit behavior for small x is D*.(t) = 1/ (aT((x)T l ~ a ). Summarizing, 



1 for small T, 

rl „. for larger. 



D* c (t)={ » r («K'-«' (15) 



a 



|r(-o)|T> 

So, two power laws coexists, with exponents I — a and 1 + a. 

Which is the basin of attraction of the fixed-point distribution, D* (t)? Let us consider D c (%) ~ 
A/t 1+cc for large T, with < a < 1. Among the class of functions that do not have a finite mean, 
these are by far the most important. We already know that the Laplace transform of D c (t) behaves 
as D c (s) = 1 +AT(-a)s a + . . ., for s — > and so, D c (p r s) = 1 - ap ra s a + .... Substituting into 



the new equation for TD c (s), and taking into account that p — > 0, 

pD c (p r s) p(l-as a p ra + - 



TD c (s) 



1 -qD c (p r s) 1 - (1 -p){\-as a p ra + ■ ■ •) 
l-asV' a + - l-as a p ra + --- 1 



+ • 



l-(l-p-as<V" a + ---) l+a5 a /j m - 1 + --- \+as a 
where we have used that r = 1/a, which, remember, means that the rescaling depends on the 
power-law exponent of the waiting-time density. Rescaling in this way ensures the existence of an 
attractor for the waiting-time densities that behave as a power law for long times, and this attractor 
is given by the fixed point (fT5l) . 

Why does this counterintuitive rescaling provide the fulfillment of a scaling law? The reason is 
in the generalized central limit theorem. Imagine that p = 1/2, so we remove half of the events; in 
order to recover the original pattern (characterized by the same waiting-time probability density) 
we need to multiply the time interval under consideration not by l/p = 2 but by 1/ p l / a , which 
is larger than l/p for < a < 1 . This is due to the fact that as the process evolves in time, 
longer waiting times appear, due to their power-law tail, being necessary to consider far longer 
time intervals. 



Reference Ib2ll puts numbers into this simple intuitive explanation. The total time interval up to 
the N—th event is fa = Ln=l T «- The largest waiting time X m (N) among the N values of the waiting 
time can be estimated as N J^ N ^D c (x)dx ~ 1, this yields x c (N) N l ^ a for large A^, and this 
means that during the time interval of lenght f# the process "does not see" the tail of D c (x) beyond 
x m (N) and one can effectively truncate the distribution at this value. Then, the "typical" value of 
fa can be associated to the mean value of the truncated distribution, so % ~ N(x) 

trunc ^ 

N l l a , and 

in this way the total time up to the N—th event scales in the same way as the largest waiting time. 

Indeed, the generalization of the central limit theorem introduces rigor in this argument, see for 
instance Ref. [j52Q. If the T n 's are power-law distributed, with a tail D c (x) =A/x l+a , the standard 
central limit theorem does not hold for t n and one has to apply its generalization. Rescaling as 
Zn = tft /N l / a the theorem states that the variable Zyy has as a limit distribution, for Af — > °°, one of 
the so called Levy stable distributions (whose mathematical form is not relevant for our purposes). 

Coming back to our problem, there is a particular value of a for which the exact inverse Laplace 
transform of the fixed-point distribution can be easily obtained. Indeed, if a = 1/2, we get that 



zr c (T) = -=-- 5 -erfc(^-) ) (16) 



whose inverse Laplace transform turns out to be 

1 e^if ^fy/r 

a\pKT a 

where erfc(x) is the complementary error function, erfc(x) = 2n~ 1 / 2 ' dx [51]. Its asymp 
totic expansion will be useful to learn how D*(t) behaves, 



e- %1 ( . 1 



QMx) = 7TxV-2xt + - 



this leads to 



C 1 



2^T 3 /2 T 3/2' 

for large T, with a = — AT(— 1/2) = 2y/nA. This is indeed coincident with the tail of the original 
distribution. On the other hand, for small argument, erfc(jc) ~ 1 and then 

D*(t) ~ l —=, 

cK ' 2An^x 

for short times, in agreement with the previous result for a general value of a, Eq. (fl5|). 
We can test these results simulating the process for instance for 

° M = nTT^ <17) 

for which A = al a . The results appear in Fig. [TJ Figure [I] (a) shows this distribution, and the 
distributions that result after rising the size threshold, keeping 10 %, 1 %, etc. up to 0.01 % of 
the events. Figure [H(b) displays the same distributions but including the rescaling, which allows 
to see the complete effect of the transformation (fT2~l) for different values of p. It is clear how 
for p = 0.01 the distribution is very close to the expected fixed-point distribution (fT6l) , and the 
agreement improves for smaller p. 

Can these results be useful for crackling-noise systems or other complex systems? Solar flares, 
when they are close to the minimum of the solar cycle 111(1 . and also e-mail activity from indi- 
viduals |53] and human-movement episodes [54] show power-law distributed waiting times, with 
exponents 1 + a about 1.5, 1, and 1.8, respectively. However, there are no indications of a second 
power law with exponent 1 — a, as requested by our theoretical calculations. 

Nevertheless, these are peculiar systems; first, in the case of solar flares, the thresholds that 
define the waiting-time distributions are not size thresholds but intensity thresholds (size is defined 
as the integral of intensity over time), and it is not clear how this change can modify the properties 



of the system. In addition, the scaling which defines the scaling law is done with the mean rate, 
as it would correspond to a distribution with a well-defined mean. Second, for human movements 
and e-mail activity, the scaling law is not achieved by means of a thinning of the process through 
the increasing of a threshold, rather, individuals with different rates are compared. And as a fourth 
example we could consider the BTW sandpile model, whose behavior (and the approach with 
which it has been studied) is very similar to that of solar flares, with a waiting-time exponent 
around 1.7 [421. 



D. Processes with Short-Range Correlations 



Reference 114911 introduced a very simple point process in which each waiting time was corre- 
lated with the size of the previous event, in such a way that waiting times following larger events 
were drawn from a waiting-time density with a short characteristic time and waiting times after 
small events had a longer characteristic time. The transformation of the waiting time density when 
the threshold is increased and the time is rescaled verifies an equation which is a generalization of 
the previous cases, 

TD c (s) 



pD c ^p r s) 



\-D c (p r s)+pD c] (pr s ) 

where D c ^(s) is the Laplace transform of D c |(t), which is the probability density of the waiting 
time that follows an event of size c' > c; more precisely D c ^(t) = D c {Ti\si-\ > c'). 

As in the last subsections, for p — > we can write D c (s) = 1 — as a + ... and D c ^(s) = 1 — 
a's a ' + . .., where a', which gives the scale of the distribution, may depend on p, increasing as 
p — > 0. Substituting in the equation for the transformation T, 

TD c (s) — P(1 -^V«' + -) 



1 - ( 1 - as a p ra + . . .) + /?( 1 - a's a 'p ra + ■■■) 
\-a!s a 'p ra ' + ••■ 1 - a's a 'p r a ' + ■ ■ ■ 1 



+ • 



y as a p ra + ---+p(l -a's a 'p ra + ---) l+as a p ra - 1 + ■■■ \+as a 
with, as usual, r = 1/a; so, whatever the dependence of a' on p, the shape of D c ^(t) is totally 

irrelevant to determine the asymptotic behavior of the process, provided that a'p — ► when p — > 0. 

This demonstrates that for the simple model we are considering, short-range correlations are not 

enough to escape from the attraction of the renewal fixed-point distributions. For the case of 

processes with a finite mean, a demonstration was already provided by Molchan, but the author 

believes the one here is more direct 



IV. DISCUSSION 



We have seen how, when a mean waiting time exists, a process without correlations cannot 
account for the appearance of a scaling law for waiting-time distributions (except in the trivial 
case of a marked Poisson process, which is not observed in real systems). In the same way, when 
a process without correlations is characterized by the absence of a finite mean waiting-time, the 
theoretical results do not compare well with observations either. Therefore, correlations build the 
shape of the waiting-time distribution in crackling-noise as well as non-crackling noise systems 
(at least for the cases studied so far, see Sec. [II]). 

But short-range correlations do not seem enough to break the dominance of the trivial Poisson 
fixed point when a mean exists, or the double power-law distributions expected when the mean 
does not exist, as we have shown for the simple example introduced in Ref. Il49n . As in equilibrium 
critical phenomena 115 611 . long-range correlations should be necessary in order to escape from the 
basin of attraction of the trivial fixed point. Further research is necessary regarding this issue, both 
from a fundamental point of view and with the goal of finding stochastic models of the systems 
displaying scaling laws. 

A promising approach is that of Lennartz et al. [57], where a long-range correlated series of 
magnitudes is generated, associating each magnitude value to a discrete time (the authors have 
in mind earthquakes, but the results are more general). Then, only extreme events, i.e., events 
above a large magnitude threshold, are considered, and the corresponding waiting-time statistics 
is obtained, with a surprising agreement with observational earthquake data [37]. In other words, 
starting with a delta distribution of waiting times for the initial process, renormalization leads to a 
nontrivial fixed point. It is an open question why this is so. 

The fact that a nontrival (nonexponential) scaling law for wai ting -time distributions may exist 
has been criticized by Molchan [58] and Saichev and Sornette [59]. Assuming that seismic oc- 
currence is well described by the ETAS model, the latter authors were able to derive the form of 
the waiting-time probability density. Here we just mention that the ETAS (epidemic type after- 
shock sequence) model is the simplest modeling of the process of earthquake triggering that puts 
all earthquakes on the same footing: any earthquake triggers other events with a probability that is 
proportional to two main factors: the Omori law, which controls the decay in time of the seismic 
rate, and the productivity law, which links the rate with the magnitude of the triggering earthquake. 
In addition, the magnitudes of the resulting triggered earthquakes are drawn independently from 



the Gutenberg-Richter distribution. Despite its simplicity, the mathematical treatment of the ETAS 
model becomes an authentic tour deforce. In any case, Saichev and Sornette get that D c (t) fulfills 
something like a scaling relation, Eq. (0Q), but with a scaling function which is not a only function 
of the rescaled waiting time, 



F{x,e) 



n6e u ( n£' 



exp —(1 — n)x- 



1-0 



where 1 + 9 is the exponent of the decay of the rate with time, given by the Omori law (6 ~ 0), 
n < 1 is the branching parameter, defining the number of events triggered directly by any event, 
and £ = AC, with C the small time constant that avoids divergence at zero time in the Omori law. 
So, in addition to the dependence on the rescaled waiting time x = X T, the density depends on the 
rate A through £, in contrast with the idea of scale invariance. 

The fact that the ETAS model does not fulfill an exact scaling law does not see m highly sur- 



6j. At 



prising, after all, due to the fact that the ETAS model itself is not fully self-similar [60, 
the end, for vanishing rate, e — > 0, it turns out that F(x) tends to an exponential distribution, and 
then the ETAS model renormalizes into the trivial Poisson process. It would be of the maximum 

nn 

interest to calculate if real self- similar models, as the Vere- Jones model H60L 16111 or the Lippiello- 
Godano-de Arcangelis model Il62ll fulfill the scaling relation (Q~J) and which is the corresponding 
scaling function. This is of course an unsolved problem. 

As a final comment, let us mention that the models used in this paper are purely stochastic, or 
mathematical. Some readers, however, could ask for a more physical approach. For the author, 
the situation is analog to the study of diffusion processes using random-walk models [5_2fl: the 
outcome is robust and independent on small details about the nature of interactions (molecular 
collisions in one case and event-event triggering mechanisms in the other). Nevertheless, there are 
examples of crackling systems that have been successfully modeled on physical grounds, using 



for instance the random field Ising model y] 



, the so called ABBM model [63 



6411 . or models of 



dislocation dynamics in plastic deformation [65]. It would be of the maximum interest to explore 
the time structure of events in those models. Regarding the geosciences, which are the systems we 
have in mind for this paper, the situation is more complicated, as the physics of those phenomena 
is still poor understood, and controlled experiments cannot be performed. It is a great challenge to 
find microscopic models of natural catastrophes that give rise to the self-organized structures that 
emerge in the long-run limit in those systems. 



The author thanks M. Boguna for his quick guide to Abelian and Tauberian theorems. 
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FIG. 1 : Illustration of the thinning plus rescaling transformation. A marked renewal process is simulated, 
with waiting-time distribution given by Eq. (fTTT ). using a = 0.5 and £ = 10. (a) The effect of thinning, which 
removes events with probability 1 — p is shown for different values of p. (b) The complete transformation T, 
adding rescaling by p 2 to thinning, shows how the resulting distribution approaches the fixed-point solution, 
Eq. Cl6]>. 



